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Chapter  1 
INTRODUCTION 

The  properties  of  metal-oxide-semiconductor  (MOS)  field 
effect  transistors,  namely,  high  circuit  packing  density,  low 
power  dissipation,  simplified  processing,  and  therefore,  low 
cost,  have  made  Very  Large  Scale  Integration  (VLSI)  a  reality. 
The  integrated  circuit  in  which  thousands  of  transistors  may 
be  fabricated  in  a  single  monolithic  chip  has  added  a  new  di¬ 
mension  to  electronic  circuit  analysis.  In  addition,  the 
complexity  and  long  term  reliability  requirements  of  VLSI 
chips,  such  as  RAM's,  require  their  design  to  be  very  nearly 
optimum  in  relation  to  both  circuit  and  device  processing. 

Also  the  designs  of  custom  logic  array  demand  a  high  degree 
of  accuracy.  One  reliable  method  of  achieving  the  above  de¬ 
sign  goal  is  through  the  use  of  circuit  simulation  programs. 
The  correct  correlation  between  the  computer  simulated  re¬ 
sults  and  experimental  measurements,  however,  depends  primar¬ 
ily  on  two  important  factors.  The  first  is  the  accuracy  of 
the  computational  method  used  in  the  simulation;  the  second 
is  the  accuracy  of  the  modeling  technique  used  in  character¬ 
izing  the  devices  in  the  circuit. 

As  the  number  of  transistors  in  the  circuit  increases, 
simulation  which  can  accurately  predict  circuit  performance 
becomes  very  expensive.  The  two  major  limitations  are  the 
computer  storage  capabilities  and  computer  time  requirements. 
Therefore,  the  determination  of  the  solutions  in  circuit 


analysis  is  important  not  only  from  a  circuit-theoretic  point 
of  view  but  even  more  so  from  a  computational  point  of  view. 
This  research  is  concerned  with  the  development  of  a  piece- 
wise-linear  model  and  some  waveform  approximation  techniques 
for  fast  and  relatively  accurate  time-domain  simulation  of 
MOS  LSI  circuits.  The  results  are  found  to  be  comparable  in 
accuracy  to  those  obtained  by  SPICE2  [iL  which  in  fact  will 
be  used  for  comparison  purposes. 

The  iterative  piecewise-linear  (PWL)  method  [21  is  used 
for  analysis  and  modeling  of  the  MOS  devices.  The  basic 
philosophy  of  the  iterative  PWL  method  consists  of  treating 
a  nonlinear  network  as  if  it  was  a  number  of  distinct  lin¬ 
ear  networks  having  identical  network  topology.  Each  linear 
network  can  then  be  analyzed  independently  by  simple  linear 
techniques.  This  is  possible  because  by  assumption  the  i-v 
curves  of  all  nonlinear  elements  can  be  realized  by  PWL  seg¬ 
ments.  Based  on  this  philosophy,  a  PWL  model  for  the  MOS 
device,  which  consists  of  an  interconnection  of  two-terminal 
elements  [3]  , is  constructed  and  studied. 

As  a  result,  the  model  as  presented  here  is  in  a  form 
most  suitable  for  easy  table  lookup  technique.  The  advantage 
of  generating  tables  of  device  characteristics  prior  to  the 
analysis  is  to  replace  the  expensive  analytic  function  eval¬ 
uations  with  simple  table  lookups.  Although  these  tables  may 
occupy  more  memory  space  than  the  analytic  functions,  this 
penalty  can  be  justified  in  terms  of  the  speedup  in  computa¬ 


tion  . 


The  technique  of  PWL  function  generation  involves  ap¬ 
proximating  the  i-v  curve  of  the  model ' s  two-terminal  ele¬ 
ments  with  continuous  straight-lines  segments.  It  is  easy  to 
see  that,  by  increasing  the  number  of  segments,  a  PWL  function 
can  represent  the  actual  curve  to  a  high  degree  of  accuracy. 
However,  more  segments  will  be  penalized  on  more  memory  space 
for  more  table  entries  and  longer  computation  time  for  more 
breakpoint  detection  [4].  For  practical  purposes,  it  has 
been  found  that  a  three-segment  PWL  approximation  gives 
acceptable  accuracy  in  digital  circuit  simulation.  Since 
there  is  no  general  method  to  find  the  best  approximation  of 
the  i-v  curves  by  segments,  the  input-output  characteristics 
of  an  inverter  will  be  used  to  investigate  the  accuracy.  The 
Katzenelson  technique  [5,  6]  will  be  employed  for  obtaining 
the  DC  operating  points.  The  transient  behavior  of  the  in¬ 
verter  is  studied  for  model  validation.  A  multiple-input 
NAND  gate  is  used  to  study  the  effects  of  a  floating  node 
caused  by  stacking  up  transistors.  Initially,  a  step  input 
waveform  is  applied  during  the  analysis  to  study  the  effects 
of  the  MOS  model  parameters  on  the  accuracy  of  the  results. 

It  is  well-known  that  the  step  response  of  a  linear  first- 
order  RC  circuit  is  exponential  with  the  time  constant  equal 
to  RC  [7]  •  This  simple  observation  is  the  key  to  the  devel¬ 
opment  of  simulation  techniques  developed  in  this  thesis. 

In  practice,  the  constant  step  input  is  not  very  realistic. 

In  most  cases,  the  signals  are  not  step  functions,  which,  as 
will  become  obvious  later,  tend  to  slow  down  the  simulation 
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method  proposed  in  this  thesis.  For  this  reason,  two  methods 
are  used  to  approximate  general  waveform  functions:  discrete 
stepwise  approximation  and  PWL  continuous  (ramp)  approxima¬ 
tion.  A  Gauss-Seidel  Waveform  Relaxation  (WR)  method  [7]  will 
be  used  to  decouple  the  circuit  equations  into  a  sequence  of 
first-order  piecewise-linear  differential  equations. 

Chapter  2  describes  briefly  the  PWL  transistor  model  and 
its  features  for  table  lookups.  Chapter  3  discusses  the  func¬ 
tion  generation  by  PWL  approximation.  Chapter  4  examines  the 
transient  behavior  of  the  model  and  presents  simple  solutions 
of  the  step  response  in  normal  form.  Chapter  5  illustrates 
the  results  with  simple  stepwise  waveform  representation  and 
suggests  an  alternative  approach  to  represent  the  waveform 
which  could  lead  to  variable  step  size  approach.  In  the 
final  chapter,  conclusions  are  presented  and  areas  for  future 


work  are  suggested. 


Chapter  2 

DEVICE  MODELING  FOR  MOS 
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2 . 1  Introduction 

In  circuit  simulation,  a  simple  MOS  model  capable  of 
producing  accurate  solutions  at  high  computational  speed  is 
very  desirable.  Hence,  the  device  physics  level  modeling 
is  impractical  for  the  application  in  VLSI  circuits,  even 
though  in  principle  the  resulting  circuit  model  would  be  ex¬ 
tremely  realistic.  Although  no  general  theory  of  device 
modeling  is  presently  available,  the  transistor  DC  character¬ 
istics  are  usually  represented  by  nonlinear  polynomial  (quad¬ 
ratic)  functions  [8].  This  representation  will  be  referred 
to  as  the  "functional"  model.  In  this  thesis,  the  character¬ 
istics  will  be  approximated  by  a  PWL  continuous  function. 

This  approximation  will  be  referred  to  as  the  PWL  model. 

In  Section  2.2,  the  "functional"  model  is  presented  which 
will  form  the  basis  for  the  development  of  the  PWL  model.  In 
Section  2.3,  the  relationships  between  the  configuration  of 
the  model  and  table  lookups  are  discussed.  In  Section  2.4, 
the  elements  of  the  PWL  model  are  shown. 

2 . 2  Functional  Model 

Without  loss  of  generality,  the  NMOS  transistor  will  be 
used  for  illustrative  purposes  throughout  this  thesis.  The 
basic  Sah  model  [8]  is  adopted  as  the  functional  model  which 


gives  analytical  expressions  for  the  NMOS.  This  simplifed 
version  of  Sah's  model  has  been  derived  from  a  circuit  de¬ 
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sign  standpoint,  considering  the  device  parameters  which 
are  important  in  the  design  of  MOS  integrated  circuits.  As 
a  practical  matter  of  VLSI  circuit  simulation,  normal  values 
and  statistical  variations  of  the  threshold  voltages,  body- 
effect  coefficient,  and  temperature  sensitivity  are  not  in¬ 
cluded.  The  special  geometrical  feature  of  a  short-channel 
device  usually  requires  a  separate  model  from  the  long-channel 
device;  this  model  does  not  include  the  short-channel  effect. 
The  device  symbol  which  represents  an  n-channel  enhancement 
type  transistor  is  shown  in  Fig.  2.1. 

The  mathematical  equations  for  the  channel  current  in 
the  triode  (non-saturation)  and  saturation  regions  are  given 
as  follows: 

Nonsaturation  Region: 

rDS  =  K  [2<VGS-VVDS-VDS2]  0  ~  VDS  VGS-VT  (2‘1) 


k(vgs-vt' 


(2-2) 


Drain 


The  constant  K  is  given  by 


K 


ue 


2t 


ox  W 
ox  L 


(2-3) 


where  p  =  average  surface  mobility  of  carriers  in  channel 
(electrons  in  n-channel)  devices 
e  =  permittivity  of  the  oxide  insulting  layer 

O  2\ 

t  =  thickness  of  oxide  under  gate 
ox 

L  =  length  of  channel 
W  *  width  of  channel 

The  VGg,  vDg»  and  V^,  are  gate  to  source  voltage,  drain  to 
source  voltage,  and  threshold  voltage,  respectively.  Two  para 
meters  are  defined  as  KP  and  k  for  convenience.  The  trans¬ 
conductance  parameter  KP  is  defined  as 

ueox 

KP  =  — - —  (2-4) 

rox 

The  width-to-length  ratio  of  channel,  k,  is  given  by 


If  the  table  lookup  method  was  used  to  solve  the  given 
mathematical  equations,  a  straightforward  approach  would 
require  a  2-dimensional  table,  Td,  to  obtain  the  channel  cur¬ 
rent  . 


V  =  V  -  V 
VGSE  GS  T 


^S^GS  ^DS'*  =  Td(VGSE,VDS) 


( 2-6a ) 
( 2-6b ) 


However,  the  number  of  entries  for  the  2-dimensional  table 


may  offset  the  advantage  of  evaluation  time.  One-dimensional 
tables  which  reduce  both  memory  requirements  and  model  evalu¬ 
ation  time  are  proposed  in  the  next  section. 


2 . 3  Equivalent  Circuit  Model  w+th  2-Terminal  Elements 

As  shown  later,  1-dimensional  tables  achieve  considerable 
memory  savings  as  compared  to  2-dimensional  tables.  One  ap¬ 
proach  is  to  use  two  1-dimensional  tables  f1-']  in  the  follow¬ 
ing  form: 

IDS  (VGS ' VGD^  -  +  Tf(VGD>  {2‘7) 

The  above  form  is  feasible  by  noting  that  the  nonlinear 
relationships  between  channel  current  and  the  terminal  volt¬ 
ages  can  be  written  as  the  difference  between  two  nonlinear 
current  components,  each  of  which  is  dependent  on  only  one 
pair  of  terminal  voltages  [lO]. 

IDS  *  J1  "  H  (2'8) 

where  I±  =  K  *  f(VQ,vs} 

I2  =  K  *  f<VG,VD) 

The  form  with  two  current  components  is  easily  generated  by 
substituting  VDS  with  VGS“VGD  ^nto  E<3*  (2-1)  .  The  final 
form  is  given  by  [10]. 

IDS  "  K(VGS-VT>2  -  K(VGD-V2  (2-9) 

Thus,  f(V_,Vx)  has  the  following  conditions: 


f(VG’VX}  =  [(VG~V  '  V  f°r  VG_VT  >  VX  (2-10a) 


for  VQ-VT  <  Vx 


(2-10b ) 


The  configuration  of  two  current  sources  is  illustrated  by 
the  model  of  Fig.  2.2. 

Proceeding  one  step  further,  an  approach  dealing  with 

interconnections  of  2-terminal  circuit  elements  yields  the 

Ebers-Moll  "like"  model  of  Fig.  2.3.  Both  forward  and  reverse 

a's  must  equal  unity  in  order  for  the  gate  current  1=0,  as 

s 

it  must  be  for  an  MOS  transistor.  The  i-v  characteristics  of 
the  diode  element  are  as  described  by  Eq.  (2-8)  and  Eq.  (2-10) 
The  given  conditions  of  f(VQ>vx)  enable  Eq.  (2-8)  to  represent 
the  channel  current  in  both  regions  with  I2  =  0  in  the 
saturation  region  and  exhibit  the  familiar  quadratic  character 
istics.  The  completely  bidirectional  feature  of  the  model 
also  reduces  the  two  1-dimensional  tables  to  one. 


2.4  Elements  of  PWL  Model 

A  one-dimensional  PWL  function  can  be  represented  math¬ 
ematically  by 


y  =  m  .  x  +  k  .  a  .  <  x  <  b  . 
J  J  J  -  -  J 


j  =  1,2, .. . 


(2-11) 


The  linear  function  defined  in  each  interval  of  Eq.  (2-11) 

will  be  called  segment,  and  the  endpoints  of  each  segment 

will  be  called  breakpoints.  A  set  of  four  parameters  (m.,k., 

J  J 

a.,b.)  is  required  to  describe  each  segment.  After  the  setup 
1  J 

of  the  PWL  segments,  an  equivalent  circuit  consisting  of  linea 
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elements  and  independent  sources  could  be  used  to  represent 
the  curves  in  a  given  region. 

The  i-v  curve  characterized  by  Eq.  (2-10)  will  be  approx¬ 
imated  by  PWL  segments.  Let  g.  be  the  slope  of  segment  j 

J 

in  the  PWL  representation.  Let  segment  j  or  its  extension 

intersect  the  i-axis  at  I..  Let  segment  j  be  defined  for  all 

a.  <  v  <  b..  Within  the  prescribed  intervals  for  segment  j, 

J  J 

the  characteristics  can  be  replaced  by  the  Norton  equivalent 
circuit  of  Fig.  2.4.  The  conductance  is  used  to  realize  the 
slope  of  the  segment  and  the  independent  current  source  is 
used  to  represent  the  i-axis  intercept.  In  order  to  keep  the 
gate  current  1^  =  0,  as  it  must  be  for  MOS,  two  linear  control 
led  current  sources  are  added  to  implement  this  physical 
property.  As  a  result,  a  PWL  model  is  constructed  as  shown 
in  Fig.  2.5.  It  is  also  noted  that  the  completely  bidirec¬ 
tional  feature  of  the  MOS  devices  gives  the  identical  i-v 
characteristics  for  the  two  nonlinear  resistors.  Hence,  the 
parameters  needed  to  describe  the  segments  are  tabulated  only 
once . 


i 
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Chapter  3 

PIECEWISE -LINEAR  FUNCTION  GENERATION 

3 . 1  Introduction 

The  linear  function  defined  In  Eq.  (2-11)  formulates 
many  segments  and  breakpoints.  It  is  easy  to  see  that,  by 
increasing  the  number  of  segments,  a  piecewise-linear  function 
can  represent  the  actual  function  to  a  high  degree  of  accura¬ 
cy.  However,  four  parameters  are  required  to  describe  each 
segment.  More  segments  will  require  more  memory  space  for 
more  table  entries  and  also  longer  computational  time  for 
more  breakpoint  detection.  For  digital  circuit  applications, 
we  found  that  a  three-segment  PWL  representation  with  appro¬ 
priately  chosen  breakpoints  gives  acceptable  solutions  when 
compared  to  SPICE2 . 

Without  any  loss  of  generality,  the  standard  depletion¬ 
load  inverter  will  be  used  for  accuracy  measurement.  The 
device  symbol  for  the  inverter  is  shown  in  Fig.  3.1.  The 
following  parameters  of  the  inverter  will  be  used  throughout 
this  thesis.  Let  kQ  and  k^  be  width-to-length  ratios  of  the 
channel  for  the  driver  and  the  load,  respectively.  Let  VTD 
and  VTL  be  the  threshold  voltage  for  the  driver  and  the  load, 

respectively.  The  transconductance  parameter  KP  is  set  to 

2 

the  default  value  of  SPICE2  at  20  /'A/V  .  Typical  values  of 
1  V  and  -2  V  are  used  for  the  threshold  voltages  of  the  en¬ 
hancement  and  depletion  transistors,  respectively.  The  power 
supply  V  ,  is  5  V.  Base  on  the  minimum  gate  area  criterion 
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and  an  arbitrary  value  of  0.2  V  for  solid  logic  "0",  kQ  and 
have  values  of  1.60  and  0.625,  respectively.  The  PWL 
model  for  the  inverter  is  shown  in  Fig.  3.2. 


L-. 


3 . 2  Approximation  by  Chord 

One  straightforward  approach  to  generate  a  piecewise- 
linear  function  is  to  connect  a  set  of  sampling  points  with 
straight  line  segments  [4].  The  i-v  curves  of  the  nonlinear 
resistors  with  3-segment  approximation  are  shown  in  Fig.  3.3. 
After  the  setup  of  the  segments,  the  four  parameters  (nij,k. , 
3j,bj)  are  tabulated  in  Tables  3.1(a,b). 

In  order  to  measure  the  accuracy  of  the  segment  approxi¬ 
mation,  the  Katzenelson  techniques  will  be  used  to  find  the 
DC  operating  points  of  the  inverter.  Results  from  the  analyt¬ 
ical  method  on  the  functional  model  will  be  used  for  compari¬ 
son.  Analysis  of  the  functional  model  gives  the  switching 
point  where  input  V^n  is  equal  to  output  VQut  at  2.25  V.  In 
other  words,  VQut  is  expected  to  have  a  high  value  when  V^n 
is  less  than  2.25  V  and  vice  versa. 

The  active  range  (-2  V  to  0)  of  the  depletion  transistor 
is  relatively  short.  Hence,  two  segments  with  equal  inter¬ 
vals  provide  a  good  approximation.  More  attention  was  placed 
on  the  selection  of  the  breakpoint  V on  the  longer  active 
range  (1  V  to  5  V)  of  the  enhancement  driver.  Selecting  the 
wrong  Vz  will  result  in  unrealistic  solutions.  This  problem 
can  be  demonstrated  by  considering  the  case  where  V^n  ^  2.25  V 
and  Vz >  2.25  V.  The  expression  of  output  V  at  the  output 


VGX (volts) 


=  FWL  segments 
=  Actual  i-v  curve 

ig.  3.3b.  PUL  approximation  by  chord  for  the 
transistor . 


Depletion 
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Table  3.1a 

Parameters  of  the  Driver  Transistor’s  i-v 
Curve  by  Chord  Approximation 


k 

m^(mhos ) 

VA) 

ak  ( V ) 

bk(y) 

1 

0 

0 

-5 

TD 

2 

KD(VVTD) 

-KD<VVTD>VTD 

VTD 

V 

z 

3 

KD(Vplus+V2VTD> 

KD<VTD2-Vplus*V 

Vz 

V 

plus 

Table  3.1b 

Parameters  of  The  Load  Transistor's 
i-v  Curve  by  Chord  Approximation 


k 

m^Cmhos ) 

jk(A) 

ak(y) 

D* 

tv 

< 

1 

0 

0 

-5 

VTL 

2 

-0 . 5Kl*  Vtl 

°-5Kl*vtl2 

VTL 

°-5vtl 

3 

-i.5kl*vtl 

2 

K  *  V 

TL 

0.5Vtl 

0 

3 


0 


terminal  is  given  by 


I,  _-I„+idd-idl+V 


gdd+gdl 


=  Vgsl*gsl+isl 


( 3-2a) 


IDI  =  Vgsd*gsd+isd 


!  3-2b) 


In  the  Katzenelson  algorithm,  gsd  and  isd  are  kept  constant 

by  the  constant  input,  while  gsl  and  isl  are  fixed  by  the  zero 

gate-to-source  potential  at  the  load.  However,  gdl  and  gdd 

will  vary  across  regions  during  the  iterative  process.  Let 

(k) 

VQut(0)  be  the  initial  guess  of  the  output.  Let  g^  and 

( k. )  •hh 

ij  correspond  to  the  parameters  of  segment  j  at  the  k 


iteration.  Let  V 


be  the  first  iterated  solution.  A 


solid  logic  0  of  0.2  V  has  been  preset  as  the  initial  guess 
at  the  output  terminal  regardless  of  the  input. 

Because  of  the  initial  guess,  the  gdd^ ^  always  start 


in  the  same  region  as  gsd. 


regardless  of  the  input.  The 


equality  of  these  two  resistors  eliminates  the  terms  that 
are  associated  with  gdd,  gsd,  idd ,  and  isd  from  Eq.  (3-1). 

The  initial  guess  also  puts  gdl  and  idl  always  in  region  1 
with  zero  values.  For  V.  <  2.25  V,  the  output  V  .  will 
vary  across  at  least  one  breakpoint  from  the  initial  guess 
to  the  final  value.  The  next  breakpoint  for  gdl  and  gdd  with 
the  rising  VQut  is  at  Vplus+VTL(3  V)and  Vin-VTD ( <  1 . 25  V) 
respectively.  Obviously,  gdd  will  switch  regions  before  gdl 
does.  In  this  case,  the  small  values  of  idl^ ^  and  gdl^  ^ 
can  be  neglected  from  Eq.  (3-1).  The  V  .  d)  is  simplified 


in  the  following  expression. 

V  (1)  isl3 

out  =  - =— — 


(3-3) 


In  the  case  where  V.  <  V__,  gdd. 


is  in  region  1 


with  the  zero  value.  Since  isl^  does  not  have  a  zero  value, 

— ,  n  "1  (1)  i.  _  ,  _  i  . .  (1)  _  _  _  _  i  ^  _ _  i  ^ 


being  zero  puts  V 


across  the  adjacent  region. 

.  7,  (1) 


In  the  case  where  2.25  V  >  V.  >  V„,  gdd. 

in  TD  J  k 

region  2.  From  Table  3.1,  Eq.  (3-3)  is  given  by 


is  m 


kL('VTL>2 


1.5625  V" 


(3-4) 


Vz-X  V 


The  condition  for  V  .  d)  to  be  able  to  travel  to  the  next 

out 

breakpoint  is  being  larger  than  V^n-VTD ( <  1 . 25  V).  Otherwise 

\T  (1)  mi  1  1  >«afn  a  1  n  i  71  +*Via  e*amn  v*fl/Y7  An  a  e*  \7  (0)  3.nd  i  t  er  a 


will  remain  in  the  same  region  as  V 


tion  will  be  terminated  with  an  incorrect  solution.  The 


above  condition  defines  the  requirement  for  as  follows: 


%  <-VTL>2 

v  <  kd 


1.5626  V*  +  1  V 


(3-5) 


After  Vz  satisfied  the  requirement,  the  accuracy  of  the 
DC  operating  point  will  be  evaluated  for  <  2.25  V.  At 
the  operating  point,  should  have  a  value  larger  than 

2.25  V.  The  final  value  of  VQut  puts  gdd  and  idd  in  region 
1.  The  expression  for  V  ,  (final)  is  given  by 


(final)  _  v 


plus 


+  isl3+idlk 


(final) 


(final) 


gsd.(flnal)  *(V.  -V__) 
_  j _ m  TD 

(final) 


(3-6) 


In  the  case  where  Vj_niVTD»  gsdj 


(final) 


becomes  gsd^ 


(final) 


which  means  the  term  associated  with  gsd^  can  be  neg¬ 


lected  in  Eq.  (3-6) .  It  is  also  reasonable  to  say  idl. 


(final) 


is  in  region  3  because  of  the  small  input.  The  resulting  ex¬ 
pression  for  V  is  free  of  Vz . 


v  (final)  _  v  =  5  V 

out  plus 


(3-7) 


In  the  case  where  gsd . ( f inal) in  region  2  (2.25  V  >  v.  > 

J  m 


„  ^  (final)  ...  (final)  (final) 

^TD  r  sno.  xqX^  ^it0  in  2T6c^xon  3  ^out 

>  4  V)  ,  the  following  expression  for  vout^"*’na^  is  obtained. 


V  (final)  =  V  .  -  2*kD* (Vz’VTD) * (Vin"VTD) 

out  plus 

kL*(VTL) 


=  5  V  -  *  (V  -1  V) 

3  Z 


(3-8) 


_  ..  ,  .  (final)  .  _  ..  (final) 

In  the  case  where  gsd.  is  m  region  2,  gdl. 

j  ,  K 

.  ...  (final)  .  .  -  ,,  ..  (final)  .  _  TTX 

and  idlk  are  m  region  2  (3  V  <  VQut  ^  4.0  V) , 

the  following  expression  for  vout^^na^  is  obtained. 


..  (final)  _  .. 
vout  plus 


2*k  (V  -V  )  *  (V .  -V  ) 
DV  Z  VTD'  V  in  TD; 

k  *V 
L  TL 


-vtl  (3-9) 


=  7  V  -  2.56* (V  -1  V) 


Equation  (3-8)  and  Eq.  (3-9)  indicate  is  depen¬ 

dent  on  the  breakpoint  V  .  Table  3.2  illustrates  the  result 
of  VQut  by  Katzenelson ' s  algorithm  with  various  values  of  V  . 

It  is  easy  to  visualize  the  accuracy  in  degrading  as  V  is 

z 

placed  farther  and  farther  from  the  input. 

On  the  segment  approximation,  the  breakpoints  give  the 
exact  value  on  the  curve.  The  error  is  increasing  away  from 
the  breakpoints  to  the  mid-point  of  each  segment.  For  the 
quadratic  function 
?  =  (x  -  vT)2 

it  is  possible  to  verify  that  the  maximum  errors  occur  at  the 
midpoint  of  each  segment  as  follows. 


Let  y  be  defined  as  Eq.  (2-11)  with  (irij  ,kj  ,a  ^  ,b^ )  . 


At 

x  =  aj  y 

*  "  «aj-VT)2 

( 3-lla) 

At 

*  =  bj  y 

-  “  (bj-VT>2 

( 3-llb) 

m,  .  V^a 

^  b  .  -a  . 

-  ■  (a.+b  )  -  2Vm 

3  3  -*■ 

( 3-llc) 

Let  error  e  be 
e  =  y  -  y 

(3-12a) 

=  (x-VT)2  -  (iHjX+kj ) 

|f  =  2(x-Vt)  -  m.  (3-13) 


Setting  de/dx  to  zero  gives 


m  . 


s  _i  +  V 

■>  +  V 


X 


( 3-14a) 


From  Eq.  (3-llc) , 


Vb3 


(3-14b) 


The  maximum  error  at  segment  2  is  as  high  as  100%  regard 

less  of  how  many  segments  are  used.  The  higher  values  of  V 

in  Table  3.2  have  the  effect  of  putting  V.  (2  V)  closer  to 

in 

maximum  error  point.  Therefore,  the  selection  of  Vz  is  con¬ 
sidered  as  crucial  in  the  approximation  method  by  the  chord. 


3.3  Approximation  bv  Tanaent 


An  alternative  way  of  approximating  functions  by  PWL  seg¬ 
ments  is  to  use  tangent  rather  than  the  chord.  The  i-v  curve 
of  a  nonlinear  resistor  is  shown  in  Fig.  3.4  with  the  new  3- 
segment  representation.  The  procedure  for  choosing  the  PWL 
segments  is  as  follows.  First,  the  voltage  axis  is  divided 
up  into  intervals.  Instead  of  connecting  the  endpoints  of 
these  intervals  to  form  the  slope  of  segments,  the  slope  is 
taken  as  the  slope  of  the  tangent  at  these  points.  The  slope 
at  each  end  point  is  represented  by  its  tangent  line  which 
defines  a  particular  segment  for  the  new  intervals.  Then  the 
intersections  of  these  tangent  lines  define  a  new  set  of 
breakpoints.  Again,  the  necessary  parameters  for  the  segments 
are  tabulated.  The  relationships  among  the  tangent  points  and 
the  breakpoints  will  be  analyzed. 

Let  va,Vb,  and  Vc  be  the  chosen  tangent  points  as  shown 
in  Fig.  3.4.  For 


i  =  K(VGX-VT) 


VT  <  VGX  <  Vend 


(3-15) 


the  slope  of  the  tangent  line  at  is 


di 

dVGX 


V 


u 


2K(W 


(3-16) 


The  slopes  of  the  tangent  line  at  va»v^»  and  vc  are  shown  in 
column  2  of  Table  3.3.  A  zero  value  of  conductance  is  used 
to  represent  the  horizontal  tangent  line  at  V  .  The  current 

a 

intercept  of  the  segment  containing  tangent  point  is 

kj  -  K(Vt2-Vu2)  (3-17) 

The  relationships  among  the  breakpoints  and  tangents  are  ob¬ 
tained  as  follows.  Let  segment  j  contain  the  tangent  point  V  , 

*  cL 

and  segment  j  +  1  contain  the  tangent  point  V^.  The  intersection 
of  these  two  tangent  lines  gives  the  new  breakpoint. 


mjX  +  kj  =  mj+lx  +  kj+l  (3-18) 

where  m.,m  ,,k.  and  k.+,  are  defined  as  Eq.  (3-16)  and  Eq. 
(3-17).  Solving  Eq.  (3-18)  gives  the  following 


XV  .  tv  . 

J  J+1 

rn.j_.-m. 

J+1  J 

k[(VT2-va2)-(VT2-Vb2)] 

2k[(Vb-VT)-(Va-VT)] 


V  +v 
a 


b 


2 


(3-19) 


The  breakpoints  in  terms  of  tangent  points  are  given  by  the 
following : 


2 


( 3-20a ) 


V1 

V2 


X 


V. +V 
D  C 


( 3-20b ) 


V 


3 


end 


( 3-20c ) 


It  is  noted  that  the  maximum  error  occurs  at  the  breakpoints 
as  shown  in  Fig.  3.4  and  their  relationships  to  the  tangent 
points  are  as  follows: 

e  s  £  (v  -V  )2 
el  4  '  b  T' 

e  =  -  (v  -V  )2 
e2  4  (  c  b; 

e„  =  K  (V  ,-V  )2 

3  v  end  c ' 

If  the  maximum  error  is  too  big, then  additional  tangent  lines 
can  be  added  at  the  breakpoint  where  the  maximum  error  occurs 
to  reduce  the  error. 

The  input-output  characteristics  of  the  inverter  will  be 
used  again  to  investigate  the  accuracy  of  the  approximation. 
The  selection  of  the  tangent  point  V.  is  necessary  to  be  V  < 
VT;  otherwise,  more  than  three  segments  will  be  established. 
The  selection  of  the  tangent  point  is  very  critical  because 
the  breakpoint  or  the  effective  threshold  voltage  is  deter¬ 
mined  by  its  value.  Although  the  maximum  error  e-^  is  also 
determined  by  V  ,  the  effect  on  the  DC  operation  is  not  sig¬ 
nificant.  However,  any  discrepancy  of  the  threshold  voltage 
in  timing  analysis  will  contribute  to  major  errors.  Hence,  it 


( 3-21a) 

( 3-21b ) 
(3-21c) 
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is  desired  to  put  Vb  as  close  as  possible  to  V^,  without  causing 
an  unacceptable  trade-off  in  the  maximum  error  eQ.  V  is 
selected  as  a  compromise  to  reduce  both  e ^  and  e^.  The  above 
criterion  of  selecting  tangent  points  is  applicable  to  both 
the  enhancement  and  depletion  transistors.  The  associated 
parameters  are  entered  in  Table  3.3  in  terms  of  the  tangent 
points . 

Table  3.3 

Associated  Parameters  for  The  Segments  by  "Tangent"  Approxi- 


mat  ion  *- 

k  m^(mhos) 

Jk(A) 

ak(V) 

Vv> 

1  0 

0 

-V 

plus 

0.5(VT+Vb) 

2  2K(Vb-VT) 

K(VT2-Vb2) 

0.5(VT+Vb) 

°.5(vc+vb) 

3  2K(V  -V  ) 

c  1 

k(vt2-vc2) 

0. 5( V  +V  ) 
c  b 

Vend 

•»  - 

Accuracy  measurements  are emphasized  on  the  approximation  of  the 
driver  because  of  its  longer  active  region.  An  arbitrary  set 
of  tangent  points  based  on  the  above  criterion  is  selected 
for  the  depletion  load.  Results  of  various  sets  of  tangent 
points  with  input  V^n  set  at  maximum  error  points  are  compared 
with  the  result  by  the  "functional"  model  in  Table  3.4.  This 
method  of  segment  approximation  provides  accurate  DC  analysis 
and  resolves  the  major  problems  of  the  chord  approximation 


method. 


Inverter  DC  Output  With  "Tangent"  Approximation 
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MODEL  VALIDATION 
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4 .1  Introduction 

Although  acceptable  DC  results  are  achieved  in  the  inver¬ 
ter  circuit,  more  detailed  study  is  needed  to  understand  the 
effect  of  the  PWL  approximation  on  the  accuracy  of  the 
transient  response.  For  an  understanding  of  the  transient 
behavior  of  the  network  model,  the  circuit  shown  in  Fig.  4.1 
will  be  analyzed.  In  addition,  the  PWL  model  will  be  used  to 
construct  a  multiple-input  gate.  The  2-input  NAND  gate  shown 
in  Fig.  4.2  is  chosen  to  study  the  effect  of  the  floating 
node  as  a  consequence  of  stacking  MOS  transistors  [ll] . 

To  simplify  the  analysis,  the  following  assumptions  will 
be  used  throughout  the  thesis: 

Assumption:  All  output  capacitances  are  combined  into 

a  simple  equivalent  load  capacitance  (CT ) .  This  sim- 

i-i 

plification  allows  an  analytical  solution  to  be  devel¬ 
oped  which  is  adequate  for  most  MOS  transient  design 
problems . 

Initially,  the  voltage  input  to  the  inverter  will  be  a 
step  function  for  both  the  turn-off  and  turn-on  cases.  The 


first  case  is  assumed  to  have  a  value  V.  =  0.2  V,  and  the 

in  ' 


second  case  to  have  a  value  of  V.  = 

m 


5  V. 
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4 . 2  Piecewise-Linear  Approach 

One  of  the  most  basic  problems  in  nonlinear  dynamic  net¬ 
work  analysis  is  to  solve  the  nonlinear  differential  equations. 
Most  of  this  tedious  process  can  be  eliminated  by  the  PWL  ap¬ 
proach.  One  way  of  analyzing  a  first-order  dynamic  network 
is  to  obtain  the  equation  of  motion  in  the  normal  form  in 
terms  of  the  state  variable  Vc  of  the  capacitor  and  solve  it. 

Because  of  the  PWL  characterization  of  the  network  model, 
it  is  possible  to  construct  at  any  time  t  =  t  a  network  which 
consists  only  of  linear  resistors  and  capacitors  since  at  any 
given  time  a  PWL  network  must  operate  on  some  linear  segment 
k.  The  response  will  continue  to  stay  in  the  same  segment 
k  until  it  reaches  one  of  the  breakpoints  of  this  segment  at 
some  time  t  =  t^.  During  the  time  interval  (t  , t^)  the  non¬ 
linear  resistance  appears  linear  to  the  energy-storage  ele¬ 
ment  (capacitor).  Therefore,  it  can  be  replaced  by  its 

Thevenin  or  Norton  equivalent  circuit.  The  resulting  network 

s  t. 

conceptually  is  the  simplest  1  order  dynamic  network  which 
has  the  general  form  of  Fig.  4.3.  The  corresponding  normal- 
form  equation  is  given  by 


c  —  +  X  =  i 

dt  R  s 


(4-1  ) 


Solving  the  differential  equation  gives  the  following  familiar 
expression , 


v(t)  =  w^vv-w^ 


—  (t  — tJ/T, 


t>tQ  (4~2a) 


■*- 


where  V  (t  )  is  the  initial  state  at  segment  k 

O 

V^(t  )  is  the  equilibrium  state  at  segment  k 

is  the  time  constant  at  segment  k 

In  the  case  where  the  PWL  segment  k  in  the  i-v  charac  »=■• 
istics  is  horizontal  or  R  +  »,  solving  Eq.  (4-1)  gives  the 
following  linear  expression: 

I 

V(t)  =  W  +  (t_to)  (4-2b ) 

If  the  above  interpretation  for  V(t)  is  applied  to  all 
times  t  >  t  ,  the  resulting  curve  will  define  the  trajectory. 
The  trajectory  always  starts  from  the  initial  state  correspond 
ing  to  t=t  .  Since  the  motion  of  the  network  is  completely 
described  by  its  trajectory  starting  from  the  prescribed 
initial  state,  the  pertinent  segment  k  containing  the  initial 
state  will  be  used  to  determine  the  equilibrium  state  and  the 
time  constant  corresponding  to  segment  k.  An  examination  of 
Eq.  (4-2)  shows  that  it  is  completely  specified  by  the  three 
parameters,  namely,  the  V,  (t  ),  V,  (t  ),  and  x,  .  Given  these 
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three  parameters,  the  solution  V(t)  can  be  obtained  by  inspec¬ 
tion  without  having  to  solve  the  associated  differential 
equation.  This  observation  is  the  key  to  the  table  lookup 
method  that  we  are  going  to  apply. 

It  is  noted  that  Eq.  (4-2a)  is  valid  only  for  all  time 

t  >  t  for  which  the  value  V(t)  falls  within  the  interval  of 
-  o 

definition  of  segment  k  and  V(t)  eventually  will  settle  down 
asymptotically  to  the  equilibrium  state  for  all  time  t  >  t  . 


If  the  value  of  V.  (t  )  is  within  the  interval  of  definition 
of  segment  k  for  all  times  t  >  t  ,  then  the  complete  solution 
consists  of  the  single  exponential  time  function.  If  the 
value  of  V,  (t  )  exceeds  the  interval  of  definition  of  segment 
k,  then  the  time  t  =  t^  where  v(t)  reaches  the  breakpoint  of 
segment  k  must  be  determined.  The  exponential  solution  is 
then  valid  only  for  t  <  t  <  t,.  From  Eq.  (4-2a),  the  time 
t^  can  be  found  easily  from  the  following  equation 


t  _  ln  w  -  vv  +  t 

1  kin  V.(t  )  -  v, )  +  lo 


(4-3) 


At  t  =  t^  V(t^)  is  chosen  as  the  new  initial  state.  The 
parameters  will  be  evaluated  at  this  breakpoint  within  the 
new  segment  and  the  process  is  repeated  as  many  times  as  neces 
sary  until  the  solution  V(t)  reaches  the  final  specified  time 
tf.  Remark:  The  DC  solution  at  t  =  0"  is  used  as  the  first 
initial  state. 


4.2.1  Inverter_circuit 

The  network  model  of  the  inverter  is  shown  in  Fig.  4.4. 
Taking  KCL  at  the  output  terminal  gives 

dV 

+  ID1  +  idl  +  gdl*vgdl  -  IL1  -  idd 

-  Vgdd*gdd  =  0  (4-4) 

where  ID^  =  Vgsd*gsd  +  isd 


IT  ,  =  Vgsl*gsl  +  isl 
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Vgdd  =  V.  -  V 
’  mo 

Vgsl  =  0 

Vgdl  =  V  -  V  . 

3  o  plus 

Rearrangement  of  Eq.  (4-4)  gives  the  Norton  equivalent  form 
of  Eq.  (4-1)  with  the  following  values  for  the  elements 

r  =  1  (4-5) 

gdd+gdl 


I  =  K,  +  K  *V. 
s  1  2  m 


(4-6) 


where  =  (idd-isd)  +  (isl-idl)  +  Vp]_us*<3dl 
K2  =  gdd  -  gsd 

In  the  turn-off  case,  the  zero  value  of  the  input  causes 
all  the  model  parameters  of  the  driver  to  be  small  which  could 
then  be  neglected  in  Eq.  (4-4) .  In  this  case,  Eq.  (4-4)  is 
governed  by  the  i-v  characteristics  of  the  depletion  load. 
Various  sets  of  tangent  points  for  the  approximation  of  the 
i-v  curve  on  the  load  have  been  used  to  observe  the  output 
risetime.  Numerical  results  using  the  PWL  approach  are  com¬ 
pared  with  those  of  SPICE2 .  Although  all  the  results  are 
close  to  those  of  SPICE2,  the  possible  accumulative  error  on 
the  delay  time  when  a  series  of  inverters  are  connected  is 
investigated.  The  set  of  (-2  V,  -1.5  v,  -0.5  V)  was  found 
to  cause  the  smallest  error  in  the  delay  time  and  thus  will 
be  used  throughout  the  remainder  of  this  thesis.  A  similar 
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procedure  is  applied  to  the  turn-on  case  to  find  a  better  set 
of  tangent  points  for  the  i-v  curve  of  the  enhancement  driver. 
Unlike  the  turn-off  case,  the  depletion  load  does  carry  some 
impact  on  the  zero-state  response  of  the  output.  Since  the 
W/L  ratio  of  the  driver  is  very  much  larger  than  the  corres¬ 
ponding  ratio  for  the  load,  the  discharging  current  through 
the  driver  will  be  much  larger  than  the  charging  current 
through  the  load.  Hence,  the  output  is  more  dependent  on  the 
discharging  current.  Initially,  when  VQut  is  high,  the  driver 
is  in  the  saturation  region  which  carries  a  large  discharging 
current.  However,  the  variation  on  the  threshold  voltage 
due  to  the  PWL  approximation  effectively  lengthens  the  satu¬ 
ration  region  of  the  driver.  As  a  consequence,  the  driver 
stays  longer  in  saturation  and  produces  more  discharging 
current  at  the  beginning  time.  As  a  result  of  the  larger 
discharging  current,  the  fall  time  by  this  model  is  shorter 
than  the  SPICE2  output.  However,  the  responses  should  be 
adequate  for  most  practical  circuit  design.  The  results  with 
the  chosen  sets  of  tangent  points  for  both  cases  are  shown 
in  Fig.  4.5. 

4.2.2  NAND  Circuit 

The  device  symbol  for  a  2-input  NAND  is  shown  in  Fig.  4.2 
and  the  PWL  model  of  the  device  is  shown  in  Fig.  4.6.  Tak¬ 
ing  KCL  at  the  output  terminal  gives 

dV 

lLi+iddi+Vgddi*gddi  =  ID^+Vgdl*gdl+idl+C^— 
where  IT .  =  Vgsl*gsl  +  isl 


(4-7) 


Tine  (ns) 

=  PWL  model  with  step  input 
=  SPICE2  output 


ig.  4.5a.  Turn-off  case  for  the  inverter 


e 


T 


T 


T 


T 


Time  (ns) 

A  =  PWL  model  with  step  input 
B  =  SPICE2  output 

Fig.  4.5b.  Turn-on  case  for  the  inverter 


IDi  =  Vgsd1*gsd1  +  isd1 


Vgddi  =  Vvo 
vgdi  =  vo-vplus 

Vgsd1  =  Va-Vf 
Vgsl  =  0 

Taking  KCL  at  the  floating  mode  gives 

ID2  +  JD3  =  vSsdi*gsdl+isd1+vgdd2*gdd2+idd2  (4-8) 

where  ID2  =  Vgdd^gdd^  +  idd-^ 

ID3  =  Vgsd2«gsd2  +  isd2 
Vgsd^  -  Va  -  Vf 

Vgdd2  =  vb-vf 
Vgddi  -  va-vo 
Vgsd2  =  0 

DC  analysis  by  Katzenelson ' s  algorithm  is  performed  by  taking 
Vj  as  an  unknown  variable  and  setting  dVQ/dt  to  zero.  Equation 
(4-7)  and  Eq.  (4-8)  are  put  in  the  following  forms  to  find 
the  DC  operating  point. 


gll 

g12 

V  I 
0 

t-i 

h-1 

- 1 

521 

g22__ 

«h' 

(4-9) 


where  g. ,  =  gdd,  +  gdl 


1^  =  isl-idl+idd^-isd^+Va (gdd^-gsd^) +Vp^us*gdl 

I2  =  isd1-idd^+idd2~isd2+Va (gsd^-gdd^) +V^ (gdd2-gsd2) 

Due  to  the  nature  of  the  logic  function  of  the  NAND  gate,  an 
initial  guess  of  5  V  is  set  to  the  output  and  the  floating 
node  for  faster  convergence.  DC  analysis  with  4  possible 
combinations  of  inputs  is  illustrated  in  Table  4.1. 

Table  4.1 

DC  Logic  Function  of  2-input  NAND 


Va 

Vb 

Vo 

Vf 

0.25  V 

5.0  V 

5.0  V 

0 . 0  V 

0.25  V 

0.25  V 

5.0  V 

2.5  V 

5.0  V 

5.0  V 

0.29  V 

0.15  V 

5.0  V 

0.25  V 

5.0  V 

3.5  V 

and 


For  the  transient  analysis, 
Vf  is  derived  from  Eq.  (4-8) 

f  gdd2+gsd^ 


the  relationship  between  V 
and  expressed  as  follows: 


o 


(4-10) 


where 


Va(gddi“gsd1  )+Vb(gsd2-gdd2  )  +  ( idd-^isd-j^ )  +  (  isd2~idd2 ) 
Vx  =  gdd^ 

Substituting  Eq.  (4-9)  into  Eq .  (4-7)  gives  the  Norton  equiv¬ 
alent  form  of  Eq.  (4-1)  with  the  following  values  for  the 
elements : 


R  = 


gsd, gdd, 

gdl+gdd^  -  gdd2+gsd^ 


=  Kl+K2*WVb 


(4-11) 


(4-12) 


where 


Kx  =  (  idd1-isd-L  )+( isl-idl  )+V plu-*gdl+  gsdl 

p  gsd-j^+gddg 

[ ( isd1-idd1 )  +  ( idd2-isd2 )  ] 
gdd2*(gdd1-gsd1) 

K2  gsd-j^+gddjg 

gsd  *(gdd9-gsd9) 

K3  ’  gsd^gdd2 

Results  with  all  the  16  possible  step  inputs  as  shown  in 
Table  4.2  are  observed.  After  the  adjustment  of  the  W/L  ratio 
on  the  load,  the  responses  of  the  NAND  gate  are  almost  identical 
to  that  of  the  inverter.  The  cases  for  rising  and  falling 
responses  are  shown  in  Figs. 4. 7  (a)  and  (b),  respectively. 


output  stays  high 

output  stays  low 

output  goes  from  high 
to  low 

output  goes  from  low 
to  high 


Voltage  (volts) 


Time  (ns) 


PWL  model  with  step  inputs 
SPICE2  output 


4.7b. 


Turn-on  case  for  the  2-input  HAND  gate 
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5. 1  Introduction 

In  general, the  signals  in  a  circuit  are  not  step  func¬ 
tions,  and  the  techniques  discussed  in  the  previous  section 
are  not  directly  applicable.  Solving  the  differential  equa¬ 
tion  of  Eq.  (4-1)  with  exponential  input,  for  example,  gives 
a  nonlinear  transcendental  equation  of  the  general  form, 

Ae_alt  +  Be“a2t  =  C  (5-1) 

Solving  this  equation  for  the  variable  t  is  a  time-consuming 
task  and  often  needs  a  number  of  iterations  to  find  the  point 
where  the  trajectory  crosses  a  breakpoint.  To  simplify  the 
task,  two  methods  of  approximating  waveforms  have  been  in¬ 
vestigated:  the  first  uses  a  discrete  stepwise  function  ap¬ 

proximation,  and  the  second  uses  a  PWL  approximation.  As 
shown  in  Chapter  4,  the  PWL  approach  with  step  input  gives  a 
simple  output  expression  which  can  be  determined  by  3  pieces 
of  information.  The  discrete  stepwise  function  defines  con¬ 
stant  level  of  voltage  at  prescribed  intervals.  During  each 
interval,  the  input  can  be  treated  as  a  step  function  which 
allows  the  analysis  by  the  method  described  in  Chapter  4.  To 
test  this  approach, a  step  input  is  applied  at  the  first  stage 
of  a  chain  of  10  inverters,  and  output  of  each  stage  is  ob¬ 
served.  The  output  waveform  of  each  stage  is  approximated  by 
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a  stepwise  function  and  fed  into  the  next  stage  as  an  input. 

In  many  cases,  when  there  is  feedback  in  a  circuit  sys¬ 
tem  the  system  will  not  be  unidirectional  and  the  subcircuits 
within  the  feedback  loop  can  no  longer  be  analyzed  as  a  se¬ 
quence  of  stages  without  iteration.  It  will  be  computation¬ 
ally  inefficient  to  consider  the  entire  circuit  as  a  big  sub¬ 
circuit.  Hence,  the  Gauss-Seidel  WR  method  will  be  employed 
to  decompose  the  circuits  into  several  dynamical  subcircuits 
each  of  which  is  analyzed  independently  for  the  entire  time 
interval  [ 7  J .  A  number  of  iterations  are  needed  for  the 
waveform  to  converge  to  the  solution,  provided  it  converges. 

A  3-inverter  ring  oscillator  will  be  used  to  illustrate  the 
application  of  WR  technique  to  circuits  with  feedback. 

Another  fact  which  may  cause  computational  problems  in 
MOS  digital  circuit  design  is  the  presence  of  pass  transis¬ 
tors.  The  pass  transistor  has  a  coupling  effect  when  the 
device  is  turned  on  and  both  its  drain  and  source  are  charg¬ 
ing  or  discharging  at  the  same  time.  A  subcircuit  which  in¬ 
cludes  pass  transistors  results  in  2  or  more  dimensional 
differential  equations.  The  WR  method  will  be  used  again 
to  decompose  the  pass  transistor  into  individual  subcircuits. 
A  simple  network  with  pass  transistor  will  be  used  as  an 
example . 

Although  the  WR  method  is  found  to  be  efficient,  it  does 


suffer  from  a  convergence  problem  when  strong  feedback  exists 
Such  a  problem  is  demonstrated  by  the  bootstrap  circuit. 


5 . 2  Discrete  Stepwise  Representation 

An  exponential  function  of  V^n  in  Eq.  (4-1)  gives  rise 
to  nonlinear  transcendental  equations.  To  avoid  solving 
these  equations,  a  stepwise  function  is  proposed  to  repre¬ 
sent  the  exponential  waveform.  A  stepwise  function  can  be 
represented  as 

Vin  =  Vj  forta<t<tb  (5-2) 

Let  V.  be  the  step  level  while  ^  and  t,  be  the  step  sizes, 
j  a  o 

During  the  interval  (td,tb),  V^n  is  constant.  Thus,  the  in¬ 
put  can  be  treated  as  a  step  function  and  solved  by  the 
method  described  in  Chapter  4 . 

A  set  of  sampling  points  is  chosen  to  set  up  the  step 
representation.  One  way  to  choose  the  sampling  points  is  to 
take  the  percentage  of  the  difference  between  the  maximum 
and  minimum  values  of  the  waveform  as  specifications  for  the 
sampling  points.  The  10-stage  inverter  chain  shown  in  Fig. 

5.1  will  be  used  as  an  example  for  the  stepwise  representation 
In  the  turn-off  case  where  the  input  waveform  is  changing 
from  high  to  low,  the  output  is  charging  up  from  low  to  high 
mainly  by  the  current  through  the  load  transistor.  A  low 
level  of  input  can  be  reached  abruptly  because  the  input 
typically  is  decaying  very  fast.  For  the  short  full  time  of 
the  input,  we  found  that  a  2-step  representation  as  shown  in 
Fig.  5.2a  gives  good  accuracy. 

In  the  turn-on  case  where  the  input  waveform  is  from  low 


to  high,  the  output  is  discharging  from  high  to  low  mainly  by 
the  current  through  the  driver  transistor.  Because  of  the 


Time  (ns) 

=  Discrete  step  function 
=  Actual  waveform 


ig.  5.2a.  Two-step  representation  of  wavefor 


short  time  response  at  the  output  due  to  large  discharge 
current  and  the  direct  impacts  on  the  discharge  current  by 
the  input,  we  found  that  a  3-step  representation  as  shown 
in  Fig.  5.2b  is  necessary  to  give  acceptable  accuracy.  Since 
the  effective  threshold  voltage  of  the  driver  is  1.5  V,  it 
was  found  the  first  sampling  point  for  step  size  should  be 
chosen  at  1.8  V  to  produce  the  turn-on  effect.  The  other 
sampling  points  will  be  determined  by  the  percentage  of  the 
difference  between  the  maximum  value  and  the  1.8  V.  If  the 
maximum  value  of  the  input  is  less  than  or  equal  to  1.8  V, 
then  the  entire  waveform  will  be  represented  with  a  step 
function  of  1.8  V.  Results  with  various  percent  points  are 
observed.  Results  of  the  inverter  chain  with  a  step  input 
and  with  the  above  representation  are  shown  in  Figs.  5.3a, 
5.3b,  5.3c  and  5.3d. 

Although  the  step  representation  allows  solving  simple 
equations,  a  considerable  amount  of  time  is  needed  to  set  up 
the  stepwise  function.  An  exponential  waveform  is  composed 
of  segments  that  are  defined  as  Eq.  (4-2)  at  a  prescribed 
time  interval.  After  a  sampling  point  for  step  size  is  de¬ 
termined,  it  is  necessary  to  locate  the  corresponding  segment 
of  the  waveform  that  contains  the  sampling  point.  After  the 
segment  is  identified,  Eq.  (4-3)  will  be  solved  to  find  the 
step  size  point.  The  setup  of  the  step  level  is  trivial  and 
straightforward.  The  other  drawback  of  the  step  representa¬ 
tion  is  the  accumulated  error  on  the  delay  time  on  the  in¬ 
verter  chain.  One  advantage,  however,  is  high-speed  computa¬ 


tion  . 


Time  (ns) 


PWL  model  with  discrete  step  waveform  representation 
SPICE2  output 


nd 


5.3a.  Output  of  2  stage. 
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=  PWL  model  with  discrete  step  waveform  representation 
=  SPICE2  output 

ig,  5.3b.  Output  of  3  stage. 


PV7L  model  with  discrete  step  waveform  representation 
SPICE2  output 

a. 

Cutout  of  6  stage. 


.  5.3c. 


?\TL  nodel  with  discrete  step  waveform  representation 
SPICE2  output 

5.3d.  Output  of  10th  stage. 


5 . 3  Network  with  Feedback 

When  a  circuit  contains  feedback,  the  subcircuits  within 
the  feedback  loop  can  no  longer  be  analyzed  as  a  sequence  of 
stages  without  iteration.  The  Gauss  -  Seidel  Waveform  Relaxa¬ 
tion  (WR)  method  [7]  is  chosen  to  decompose  the  circuit  into 
several  subcircuits  each  of  which  is  analyzed  independently 
as  an  individual  stage  for  the  entire  time  interval.  A  3- 
inverter  ring  oscillator  shown  in  Fig.  5.4  is  used  as  an 
illustrative  example. 

In  the  oscillator  circuit,  instead  of  solving  the  ordinary 
differential  equation  of  the  form 

X  =  f (X, t)  (5-3) 

as  in  inverter  chain,  a  first-order  2-dimensional  differential 
equation  with  Xf  as  the  feedback  term  in  the  form 

X  =  f (X, Xf , t )  (5-4) 

will  be  considered.  The  basic  idea  of  the  WR  method  is  to 
fix  all  the  unknown  variables  and  solve  for  one  variable  at 
a  time  as  a  1-dimensional  differential  equation.  The  solution 
obtained  for  that  variable  can  then  be  substituted  into  an¬ 
other  equation  to  solve  another  variable.  The  procedure  is 
repeated  until  convergence  to  a  solution  is  reached. 

The  step  representation  proposed  in  the  last  section  is 
not  practical  for  application  in  the  WR  method.  First,  in 
the  case  where  the  circuit  is  oscillatory,  a  lot  of  time  is 
needed  to  set  up  the  steps  for  each  rising  and  falling  edge. 
Second,  the  step  representation  also  seems  to  provide 
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insufficient  accuracy  for  nonunidirectional  circuits.  A 
step  representation  with  fixed  small  step  size  is  proposed. 
This  approach  eliminates  the  computation  of  step  size  which 
can  be  quite  time-consuming  and  the  discrete  levels  of  step 
size  still  allow  the  application  of  Eq.  (4-2) .  The  small 
step  size  also  provides  a  higher  degree  of  accuracy. 


A  fixed  step  size  of  1  nsec  is  used.  One  approach  to  set 
up  the  step  level  is  to  take  the  average  of  the  voltage 
levels  at  both  ends  of  a  step  size.  To  simplify  the  analysis 
further,  the  procedure  to  detect  the  time  point  where  the 
trajectory  crosses  a  breakpoint  is  eliminated.  Only  the 
initial  state,  the  equilibrium  state,  and  time  constant  of  Eq 
(4-2)  are  evaluated  at  the  endpoints.  Because  of  the  small 
step  size,  the  error  due  to  this  simplification  is  not  sig¬ 
nificant.  A  solution  similar  to  Eq.  (4-2)  will  be  reproduced 
in  terms  of  the  stepwise  input.  In  general,  the  first-order 


differential  equation  of  Eq.  (4-1)  can  be  expressed  as 


x  =  ax  +  bu  +  c 


where  x  is  the  output  and  u  is  the  input. 


(5-5) 


Solution  of  the  above  equation  is  obtained  by  solving 

x(t)  =  x(tQ)ea^t  to)  +  eat  j  e  aT (bu (t) +c) d^  (5-6) 

fco 

Let  the  input  u  be  expressed  as  the  average  voltage  between 
the  endpoints  of  a  step  size  as  shown  in  Fig.  5.5. 
u  (t  )  +  u  (t) 


Then 
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:-l  l 


I 


bu(t  )+c 

x(t1)  =  |  x  ( tQ )  +  - 2 -  *  h  ;  e 


.  bu(t,)+c 

+  h  (5-10) 


If  the  difference  between  u(tQ)  and  u(t^)  is  small,  the  step 
size  h  can  be  increased  without  degrading  the  accuracy.  On 
the  other  hand,  the  step  size  can  be  decreased  to  improve 
accuracy.  Results  from  the  ramp  representation  with  fixed 
step  size  shown  in  Fig.  5.8  are  almost  identical  to  those  in 
Fig.  5.6. 


The  technique  applied  to  check  for  convergence  compares 
the  values  of  the  iterated  solutions  at  the  end  of  the  given 
time  interval.  The  results  indicate  this  method  of  checking 
for  convergence  is  quite  efficient  and  reliable. 


5.4  Network  with  Pass  Transistor 

The  circuit  shown  in  Fig.  5.9  is  a  basic  part  of  the  com¬ 
mon  dynamic  shift  register.  The  PWL  model  of  the  circuit  is 
shown  in  Fig.  5.10.  This  circuit  has  a  coupling  effect  when 
the  clock  control  VCLK  of  the  pass  transistor  is  on  and  both 
and  Cj  are  charging  and  discharging  at  the  same  time. 

The  circuit  can  be  described  by  a  first-order  2-dimensional 
differential  equation  of  the  following  form: 


=  MlVCl 


M2VC. 


+  M  *V 
3  CLK 


+  VVin  + 


M, 


(5-lla) 


*  M5VC, 


M6VC- 


+  M  *V 
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+  M, 


( 5-llb) 


The  Gauss-Seidel  WR  method  will  be  used  again  to  solve  these 
equations  as  two  1-dimensional  differential  equations. 
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5.8  Output  of  the  1  inverter  m  the  three 
inverter  ring  oscillator. 


Taking  KCL  at  the  output  terminal  of  the  inverter  gives 


the  following  expression. 
dV, 


Cl 


Cl 


+  +  IT^  +  Vgdl*idl  +  idl  -  IL^  -  Vgdd*gdd 


dt 


-  idd  -  Vgdt*gdt  -  idt  =  0 


(5-12) 
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Rearrangement  of  Eq.  (5-12)  gives  the  Norton  equivalent  form 
of  Eq.  (4-2)  with  the  following  values  for  the  elements: 
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R  = 


(5-13) 


gdl+gdd+gdt 


*s  -  K1  +  VVi„  +  K3*VC2  +  K4*VCLK 


(5-14) 


where  =  (isl-idl)  +  (idd-isd)  +  (idt-ist)  + 
K2  =  gdd  -  gsd 
K3  =  gst 
K4  =  gdt  -  gst 


Substituting  Eq.  (5-13)  and  Eq.  (5-14)  into  Eq.  (5-12)  gives 
an  equivalent  form  of  Eq.  (5-lla)  which  is  also  used  as  the 
basic  circuit  equation  for  the  WR  method.  Taking  another  KCL 
equation  at  Vc2  gives 


dV 


C2 


C2 


dt 


+  Vgst*gst  +  ist  -  IT2  =  0 


(5-15) 


where  IT2  =  Vgdt*gdt  +  idt 
Vgst  =  VCLK  -  Vc2 

V^dt  =  VCLK  “  VC1 


Rearrangement  of  Eq.  (5—15)  gives  the  Norton  equivalent  form 

of  Eq.  (4-2)  with  the  following  values  for  the  elements: 

R  =  1  (5-16) 

gst 


rs  -  K5  +  VVC1  +  VVCLK 


(5-17) 


where  Kg  =  ist  -  idt 
Kg  =  gdt 

=  gst  -  gdt 


Substituting  Eq.  (5-16)  and  Eq.  (5-17)  into  Eq.  (5-15)  gives 
an  equivalent  form  of  Eq.  (5-llb)  which  is  the  basic  circuit 
equation  for  the  WR  method. 


The  delay  through  such  a  circuit  is  different  from  the 
inverter  delay  which  is  more  dependent  on  the  input  level. 
In  the  case  where  both  capacitors  are  charging  up,  the  re¬ 


sults  by  the  WR  method  show  a  slower  time  response  than  those 
of  SPICE2.  The  major  reason  is  the  larger  effective  thresh¬ 
old  voltage  of  the  pass  transistor  gives  lower  steady-state 
values  for  Vc2 .  Since  is  coupled  to  Vc2,  both  VC1  and 
Vc2  have  significantly  longer  rise  times  than  the  SPICE2  output 
The  results  are  shown  in  Fig.  5.11. 

5 . 5  Convergence  Problem  of  Waveform  Relaxation  Method 

The  bootstrap  circuit  in  Fig.  5.12  has  a  floating  capaci¬ 
tor  as  the  feedback  element.  The  WR  method  and  stepwise 
function  are  applied  to  analyze  this  circuit.  If  the  value 
of  is  relatively  large,  then  the  solution  does  not  con¬ 
verge.  A  much  smaller  step  size  is  necessary  to  analyze  this 
problem.  The  simulated  results  with  =  0.005  pF  and  the 
other  output  capacitors  at  0.1  pF  are  shown  in  Figs.  5.13a,b, 
c,  where  it  can  be  seen  that  the  result  by  this  approach  is 
not  accurate  in  comparison  with  SPICE2 . 
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5.11.  Output  of  inverter  at  the  pass  transistor 
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Chapter  6 
CONCLUSION 

A  PWL  model  of  MOS  transistor  circuits  has  been  con¬ 
structed.  The  aim  of  this  model  is  to  provide  acceptable 
accuracy  without  solving  complex  nonlinear  equations.  A 
configuration  of  interconnection  of  2-terminal  nonlinear 
elements  for  the  MOS  is  used.  This  form  not  only  satisfies 
the  requirement  of  the  PWL  canonic  form  but  also  allows  the 
application  of  1-dimensional  table  lookups.  In  addition, 
more  memory  saving  is  achieved  because  the  completely  bi¬ 
directional  property  of  the  model  requires  only  one  table  to 
represent  both  nonlinear  elements. 

Two  approaches  to  function  generation  by  PWL  segments 
have  been  studied.  Both  approaches  are  convenient  and  easy 
to  implement.  For  accuracy  measurements  of  the  approxima¬ 
tion,  the  Katzenelson  algorithm  is  used  to  study  the  DC 
characteristics  of  the  depletion-load  inverter.  The  first 
approach  uses  chord  approximation.  The  advantage  of  this 
approach  is  the  exact  representation  of  the  i-v  curve  at  the 
threshold  voltage  and  the  straightforward  setup  of  the  chords 
However,  it  suffers  from  accuracy  problems  if  the  DC  operat¬ 
ing  point  is  placed  near  the  maximum  error  point.  The  DC 
output  also  exhibits  high  dependence  on  the  endpoint  of  the 
chords.  The  second  approach  is  the  tangent  approximation. 

It  is  slightly  more  complicated  to  set  up;  however,  the  error 
has  been  reduced  significantly.  Even  with  the  DC  operating 
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point  at  the  maximum  error  point,  the  results  still  provide 
good  accuracy.  The  major  drawback  to  this  approach  is  the 
variation  of  the  threshold  voltage  sometimes  imposes  a  prob¬ 
lem  on  timing  analysis. 

For  the  understanding  of  the  model's  transient  behavior, 
step  inputs  were  applied  to  the  inverter  and  2-input  NAND 
gate  for  observation.  The  values  of  the  step  function  are 
selected  such  that  the  impacts  of  the  load  and  driver  on  the 
output  can  be  separated.  Thus,  each  transistor's  model  para¬ 
meters  can  be  examined  to  determine  the  best  PWL  segment  ap¬ 
proximation.  Although  the  changes  of  threshold  voltages  due 
to  tangent  approximation  contribute  to  errors  in  the  tran¬ 
sient  analysis,  the  results  are  considered  acceptable.  The 
solutions  of  the  PWL  model  in  the  2-input  NAND  gate  circuit 
with  various  inputs  display  accuracy  that  is  as  good  as  the 
solution  in  the  inverter. 

The  PWL  approach  allows  a  PWL  model  network  to  be  trans¬ 
formed  into  a  first-order  dynamic  network  in  a  given  segment. 
The  solution  of  this  simple  equivalent  circuit  with  constant 
input  leads  to  an  exponential  waveform  or  linear  waveform, 
whose  exact  form  can  be  found  by  determining  only  three  pieces 
of  information,  namely,  the  initial  state,  the  equilibrium 
state,  and  the  time  constant  associated  with  the  segment. 

This  information  is  easily  obtained  by  employing  the  table 
lookup  method. 

For  realistic  circuit  analysis,  the  exponential  waveform 
is  considered.  A  discrete  stepwise  function  is  proposed  to 
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represent  the  exponential  waveform.  First,  this  representa¬ 
tion  avoids  solving  the  nonlinear  equations.  Second,  it 
adopts  the  simple  solution  for  constant  input  which  can  be 
solved  easily  with  table  lookup.  A  10-stage  inverter  chain 
is  used  to  illustrate  the  effectiveness  of  the  stepwise  re¬ 
presentation.  The  outputs  at  various  stages  are  observed  to 
give  a  good  approximation  quickly.  Although  the  discrete 
representation  enjoys  solving  simple  equations,  a  consider¬ 
able  amount  of  time  is  needed  to  set  up  the  stepwise  function 
The  errors  associated  with  delay  time  are  considered  large. 

The  effect  of  feedback  in  a  circuit  is  studied  with  the 
application  of  the  WR  method.  Two  representations  of  wave¬ 
forms  are  proposed  to  be  incorporated  into  the  WR  method.  The 
first  waveform  representation  method  is  a  stepwise  function 
with  small  fixed  step  size.  The  analysis  of  this  representa¬ 
tion  is  treated  as  constant  input  at  prescribed  intervals. 

The  small  fixed  step  size  eliminates  the  computation  of  step 
size  and  gives  very  good  accuracy.  However,  the  small  step 
size  requires  a  lot  of  table  lookup.  The  second  representa¬ 
tion  method  is  a  ramp  function  which  allows  the  application 
of  variable  step  size.  It  can  result  in  less  table  lookup 
or  more  accuracy  as  described.  The  analysis  of  this  function 
is  not  more  difficult  than  for  the  step  case  because  both 
solutions  are  expressed  in  the  same  form  and  can  be  solved 
easily  with  table  lookup. 

The  WR  method  with  small  step  size  stepwise  representa¬ 
tion  was  also  applied  to  solve  the  network  with  a  pass 


transistor.  The  result  exhibits  significantly  slower  re¬ 
sponse  on  the  turn-off  case.  The  major  reason  is  the  larger 
effective  threshold  voltage  of  the  pass  transistor  gives 
lower  steady-state  value.  The  study  of  the  bootstrap  circuit 
indicates  the  WR  method  will  not  be  convergent  when  a  strong 
feedback  element  exists. 

The  simulation  data  of  the  following  cases  which  were 
analyzed  in  the  previous  sections  with  the  waveform  relaxa¬ 
tion  and  small  fixed  step  size  stepwise  representation  are 
compared  with  those  obtained  from  SPICE2. 

Case  6.1:  The  3-inverter  ring  oscillator  in  Fig.  5.4. 

Case  6.2:  The  pass  transistor  network  in  Fig.  5.9. 

Case  6.3:  The  bootstrap  circuit  in  Fig.  5.12. 

The  other  example  used  the  variable  level  discrete  steps  for 
waveform  representation. 

Case  6.4:  The  ten-stage  inverter  chain  in  Fig.  5.1. 

These  simulation  data  shown  in  Table  6.1  indicate  the 
PWL  approach  provides  some  savings  in  CPU  time.  However,  some 
examples  also  illustrate  accuracy  or  convergence  problems.  A 
modified  version  of  the  WR  method  or  a  total  new  approach 
should  be  investigated  to  overcome  the  problems.  A  different 
set  of  PWL  segments  may  be  necessary  to  represent  a  pass 
transistor.  The  waveform  representation  by  ramp  functions 
should  be  developed  in  more  detail  to  implement  variable  step 
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